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Abstract 

Open heavy flavour production at e+e^ cofliders in deeply inelastic 57 scattering has 
an interesting feature: the structure function F2{x, Q^) for this process is calculable for 
X > 0.01, and is essentially proportional to the gluon density in the photon for smaller 
X values. We give estimates for event rates at LEP2 and a Next Linear Collider in x, 
bins, and present differential distributions in the transverse momentum and rapidity 
of the heavy quark for the case of charm. We include all next-to-leading-order QCD 
corrections, and find theoretical uncertainties are well under control. 
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1. Introduction 



A copious amount of heavy quarks produced in two-photon colHsions will be generated 
at high energy e^e~ colliders such as LEP2 and a future Next Linear Collider (NLC). The 
largest fraction of these will come from so-called no-tag events, in which neither of the 
leptons generating the photons is seen. However, due to the high energy and luminosity 
of these machines, a sizeable number of them will be produced in events where either the 
outgoing e"^ or e~ is tagged ("single-tag events"). In such events, the photon coming 
from the tagged lepton is far off-shell and spacelike, hence this reaction amounts to deeply 
inelastic electroproduction of heavy quarks on a photon target. An important aspect of 
on-shell photons initiating a hard scattering is that they may behave either as pointlike 
particles ("direct" process), or as hadronic (vector) particles ("resolved" process) In 
no-tag events, both the direct and resolved process contribute to the cross section. At LEP2 
they contribute in about equal amounts to the cross section for charm production in two- 
photon collisions |^]. In such events one can separate these components by making use of 
forward detectors and using the presence of a remnant jet of spectator partons in resolved 
processes as a separator [^. The interesting feature of single-tag heavy-quark events is 
that this separation occurs quite naturally in the deeply inelastic Bjorken scaling variable 
X, as we will show below. Moreover, the presence of the heavy quark mass ensures that 
this separation is unambiguous to next-to-leading order (NLO) in QCD. The direct process 
is directly calculable in QCD and is free from such phenomenological inputs as parton 
densities in the photon. The resolved channel on the other hand is directly proportional 
to the poorly known gluon density in the photon. Therefore a reasonable sample of single- 
tag heavy-quark events allows one to confront simultaneously a well-controlled perturbative 
QCD calculation with experiment, and constrain the gluon density in the photon. 

In the past (open) heavy quark (mainly charm) production in two-photon collisions at 
e^e~ colliders has been difficult to observe in experiments due to the low charm acceptance. 
The reaction e~^e~ — > e~^e~D*^X has been thoroughly studied experimentally. The exis- 
tence of the D*^ has been inferred either from direct reconstruction |QJ5|,^] or from unfolding 
the distribution of soft pions |^J8| resulting from its decay. In addition, measurements have 
been made using soft leptons P,p^ and kaons ^ to tag the charm. JADE and TPC/Two- 
Gamma have performed experimental studies of the reaction e~^e~ e~^e~D*^X with 
one outgoing lepton tagged at low average value (Q^) of the momentum transfer squared 
of the tagged lepton (below 1 (GeV/c)^). TOPAZ [Q] has performed a study as well at 
somewhat larger (Q^). The total number of events obtained was however very small (about 
30 for TOPAZ). See ||ll| for a more extensive review of the experimental situation. 



The cross section and single particle distributions for 77 — > cc have been calculated to 
next-to-leading order (NLO) in QCD in Q, and agree with experimental results. Corre- 
lations in the direct channel were studied in |12]. For the single-tag case, the structure 
functions F2 and F2 for charm production were calculated to NLO in QCD in |13]. In this 
letter we employ the results of to estimate the expected number of single-tag events 
m X, 

Q2 bins for LEP2 and a NLC, and examine differential distributions in the transverse 
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momentum and rapidity of the charm quark. 

In section 2 we describe the formalism and explain our notation. In section 3 we show 
the structure function for charm production, give estimates of event rates in x, bins 
for the case of LEP2 and a NLC and show single-charm-quark differential distributions. We 
summarize and conclude in section 4. 



2. Formalism 



We consider the reaction (see Fig.l) 

e-(j>e) + e-(p'J + e+ + Q(pi) +X, (1) 

where Q{pi) is a heavy quark with momentum pi and X denotes any hadronic state allowed 
by quantum-number conservation laws. When the outgoing electron is tagged, this reaction 




Figure 1: Charm production in a single-tag two-photon event. 



is dominated by the subprocess 



Q{pi) + x, 



(2) 



where one of the photons is highly virtual and the other is almost on-mass-shell and trans- 
versely polarized. The case where the positron is tagged is completely equivalent. This 
process is described by the differential cross section 

mi 



dxdQ^dTidUi 



/• S 



(3) 



y 



dTidUi 



where the (fF^{x,Q'^,m?)/dTidUi {k = 2,L) denote the (doubly differential) deeply in- 
elastic photon structure functions, a = e^/47r is the fine structure constant, and m the 
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heavy quark mass. S is the center of mass energy squared of the e^e system. The Bjorken 
scahng variables x and y are defined by 

k-q , 

X = 7T1 , y = , q = Pe-Pe- (4) 

Further, we define the variables 

Ti=T -m"^ = {k-pif -m"^ , Ui = U - = {q - pif - . (5) 

The momenta of the off-shell photon and the on-shell photon obey the relations q^ = —Q^ < 
and A;^ ~ respectively. Because the photon with momentum k is almost on-mass-shell, 
eq. (3) is written in the Weizsacker- Williams approximation: the function fy{z, S/m1) rep- 
resents the probability of finding a photon 'y{k) in the positron, with longitudinal momentum 
fraction z. It is given by 

^" I, ^ '^min '*'min '*'max ) 

where k^^^ = {z^ml)/{l-z) and /cLx = (l--^) (^fe^max) V(z2m2)/(l-z). Here = VS/2 
is the lepton beam energy and ^max is the anti-tag [| angle. In the rest of the calculation 

In eq. (3) the differential structure functions can be represented as 

\ Sr^ f^^'^^ dz ^ X 2\ ir-i 



dF2{x,Q\m^) = dCk,,{x,%)+ r^-fn-,f^})dCkA^,%,rn^' 

Jx z z ' fij 



^=q,q,9 



= dF^ V dF,^'^^^ , (7) 

where d represents d^/dTidUi and Zmax = Q'^/iQ'^ + 4m?). The two terms represent the 
contributions to the structure functions from pointlike (or "direct") and hadronic (or "re- 
solved") photons, respectively. The are photonic parton densities and the dCk,i{i = 
q,q,g,^) are (differential) Wilson coefficient functions, and /i/ is the mass factorization 
scale. At lowest order, only the photonic gluon density appears in eq. (|^). F2 and F2 were 



calculated in |13] to NLO in QCD by computing all 0{as) corrections to the coefficient 
functions C^^i- We use the results of this calculation for the present analysis. Coupling con- 
stant renormalization was performed in the MS scheme, modified such that heavy flavours 
decouple in loops when small momenta flow into the fermion loop, and mass renormalization 



in the on-shell scheme. Mass factorization was also performed in the MS scheme. See [|13| 
for further details. Note the absence of the scale A*/ in the second term on the right hand 
side of eq. This is due to the fact that, through NLO, the heavy quark mass prevents 
a collinear singularity from occurring in the pointlike piece. For massless quarks this sin- 
gularity is subtracted at scale fif and absorbed into the hadronic piece. As a consequence 
F^' is calculable - the only theoretical uncertainties stemming from ag and the heavy 
quark mass- whereas F^'^^^ is primarily sensitive to the gluon density in the photon, in 

^The maximum angle below wliicli no energy deposit in forward calorimeters may be seen, tiius selecting 
events in which the lepton generating the target photon goes down the beampipe. 
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analogy to the proton case. We will see that the two components on the right hand side 
in eq. dominate in different regions of the Bjorken variable x. /,From eq. (3) we will 
estimate the numbers of events expected per x, Q'^ bin at LEP2 and a future NLC. 

In analyzing charm production, we will also show differential distributions in transverse 
momentum and rapidity of the detected charm quark in the 7*6"*" cm. frame. These 
quantities are derived from the variables defined in the above as follows. We define W to 
be the invariant mass squared 

W = {q + kf (8) 

and W = W + Q^. The transverse mass of the heavy quark in the 7*7 frame is determined 
by 

W''^m^ = W'TiUi + Q^rf + Q'^W'Ti . (9) 

The transverse mass is the same in the 7*6"*" frame. Further, in the 7*7 frame the energy 
of the heavy quark is 

2^ ^ ' 



Its longitudinal momentum can be found from the condition = E'^ — m^, so its rapidity 
in the 7*7 frame is 

In the 7*e^ frame, y = y + ^ln(z) (we define a particle to have positive rapidity when it 
travels into the hemisphere from which the 7* originates). 



3. Results 

We first list our choices for various quantities used in the calculations of this section. 
Deviations from them will be explicitly indicated. We used for m, the heavy quark mass, 
the value 1.5 GeV for charm and 4.75 GeV for bottom. We assumed a center of mass energy 

= 180 (500) GeV for LEP2 (NLC). For the strong coupling we took a one-loop running 
as with three (four) light flavours and A = 232 (200) MeV at LO and a two-loop running 
as with A = 248 (200) MeV at NLO for charm (bottom). We put the renormalization scale 
equal to the factorization Hf and chose = Q. We also imposed the following requirements 
on the angle and energy of the tagged lepton: ^tag > 30 (40) mrad for LEP2 (NLC) and 
E'tag > -E'beam/2, where -Ebeam = ^5/2. These cuts imply a minimum of 3.65 GeV^ for 
LEP2 and 50 GeV^ at our hypothetical NLC. Furthermore we used the GRV ||T^ parton 
densities for the hadronic photon results, and the Weizsacker- Williams density of ||l^ with 
a ^max = ^tag for the equivalent target photon spectrum. For the total integrated luminosity 
of LEP2 (NLC) we used / Cdt = 500 pb"^ (20 fb"^). At a NLC, beamstrahlung effects are 
expected to play an important role. We include them here by adopting for its spectrum the 
expression given in |l^, with parameters Te// = 0.039 and cj^ = 0.5 mm [^], corresponding 
to the TESLA design. When showing differential distributions we will show them from the 
point of view of the 7*e^ cm frame, i.e. we assume here that the electron is tagged and 
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Figure 2: Hadronic (dashed lines) and pointlike (solid lines) components of F2/a vs. x at 
= 10(GeV/c)^. The lower solid line is the LO, the upper one the NLO case. The lower 
dashed line at x = 0.001 is the NLO, the upper one the LO case. 



take the case where the positron is tagged into account by multiplying the results with a 
factor 2. 

We begin by showing in Fig. 2 the quantities -P^''^^ and F^'^^^ versus x both at LO 
and NLO, for = 10 (GeV/c)^. The 0{as) corrections to both quantities are for the most 
part fairly small in the x-region shown. We may therefore assume that even higher order 
contributions are negligible, and these components are thus calculated with some accuracy. 

The prominent feature in these figures is clearly the separation in x of the components. 
This fact was already noted in [13|, and we emphasize it here. As a result, the opportunity 
to confront a precise calculation and measure the small-x photonic gluon density presents 
itself in one experiment. As mentioned earlier, the difficulty of efficient charm tagging makes 
such an experimental study very difficult in practice. To judge the feasibility of such a study 
we have integrated eq. (3) for various x, bins, and obtained estimates for the number of 
charm quarks per bin produced at LEP2 and NLC by both hadronic and pointlike photons. 
Note that these estimates use the NLO calculation. 

To perform the integrals over x and Q^, we used fitted versions of the integrated NLO 
coefficient functions Ck,i{i = q,q,g,j) in eq. as the actual expressions in [13| are too 
long for fast evaluation. By adapting the fitted coefficient functions of electroproduction of 
heavy quarks on a proton target [^] to the photon target, we were able to speed up our 
computer program by as much as a factor of fifty. 

The results for LEP2 are shown in Tables 1 and 2. We give in these tables the expected 
number of events due to pointlike (Table 1) and hadronic (Table 2) photons. Assuming a 



5 



charm tagging efficiency of 1-2%, we tliat on average a few tens of events per bin sliould be 
observable for larger x values. Furthermore, to test the stability of the results we varied the 
renormalization/factorization scale fj,f from (5/2 to 2Q. The results change relatively little 
under these variations, as seen in Tables 1 and 2. For these LEP2 conditions, we varied 
in addition the charm quark mass from 1.3 GeV to 1.7 GeV and found that the pointlike 
contribution (Table 1) changed under this variation in the charm mass by less than 10% 
in all bins, except in the x bins 0.32 < x < 1. In these bins, the mass variation causes a 
variation in the position of the charm quark production threshold, and from the shape of 
F2'^^ in Fig. 2 one may readily understand that this causes a large change in the number of 
events, as much as 70% in the lowest bin. The sensitivity of the hadronic contribution 
(Table 2) to these changes is much more uniform in x, varying from 30-40% in the lowest 
Q"^ bin to 10-20% in the higher bjns. 

In Table 3 we give similar estimates for the NLC, both for charm and bottom produc- 
tion. Here we have summed the contributions due to pointlike and hadronic photons for 
each entry. We do however show separately the number of events due to beamstrahlung 
and Weizsacker- Williams bremsstrahlung contributions. Note that the number of charmed 
events here is truly large, even at large Q^. The number of events containing bottom quarks 
is much smaller due to charge and phase space suppression. 

Clearly the production rate of charm quarks produced by hadronic photons is too small 
at LEP2 to be useful for a good measurement of the gluon density in the photon. At a NLC 
the rates are significantly higher however. We therefore give in Table 4 the production rate 
for a NLC in x, bins for various choices of parton densities, charm quark mass values 
and mass factorization scales. Examining the composition of these results, we find very 
little contribution due to the hadronic channel from beamstrahlung photons (about 1-2%). 
We further find the contributions due to pointlike photons both of WW and beamstrahlung 
origin to be similarly small. There is therefore a good possibility of measuring the gluon 
density in the photon at a NLC, poor tagging efficiency notwithstanding. The most signifi- 
cant uncertainty is related to the charm quark mass, however a differentiation of the gluon 
densities of GRV ||l5[ and ACFGP |16| seems certainly feasible. 

While bottom quark production lends itself to considerably better tagging efficiencies, 
the production rate is severely suppressed relative to charm. An analysis similar to the one 
for charm reveals that a measurement of the photonic gluon density from bottom production 
in 67 scattering does not seem possible at a NLC. 

We turn to single particle inclusive differential distributions in transverse momentum 
{ptj and rapidity (y) of the heavy quark in the 7*6"*" cm. frame. Due to the low rate 
of bottom quark production we only show these distributions for the case of charm. We 
integrate x and over the intervals 0.001 < x < 1 and 3.2 < Q'^ < 320 (3200) GeV^ 
for LEP2 (NLC). To approximate detector limitations, we use the restrictions pt > 0.5 
GeV and \y\ < 2, in addition to the cuts listed in the beginning of this section. The pt 
restriction is of course quite loose, but it allows us to study the behavior of the pt spectrum 
at low Pt- Furthermore, for the results shown in Figs. 3 and 4 we took the factorization 

■)2 I ^2 



(=renormalization) scale Hf = J Q"^ + pf- 
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Figure 3: Hadronic and pointlike distributions in pt and y of the charm quark in LO and 
NLO at LEP2. a) The upper set of histograms represents the pointlike pt- distribution, the 
lower the hadronic. The solid histograms represent the LO case, the dashed the NLO. b) 
Same as in a) but now for the y- distribution. 

In Fig. 3a we present for LEP2 the LO and NLO p^-distributions of the charm quark 
for the pointhke and hadronic piece separately. In neither instance do the LO and NLO 
curves differ significantly. While the hadronic piece at NLO is larger than the LO at large 
Pt (which is similar to what was found for a proton target in ||2l|), the reverse is true for 
the pointlike piece. The bulk of the events has pt < 5 GeV. In Fig. 3b the same is shown 
but now for the rapidity distribution. Clearly the charm quark has predominantly negative 
rapidity. 

In Figs. 4a and 4b we present for the NLC the pt spectrum results due to pointlike 
and hadronic photons for a Weizsacker- Williams spectrum and a beamstrahlung photon 
spectrum. Both spectra lead to similar plots as for LEP2 in Fig. 3a. The dominance near 
small Pt of the beamstrahlung piece is due to the small final-state invariant mass preferred 
by the beamstrahlung spectrum. 

4. Conclusions 

We have argued in this letter that a measurement of the cross section for heavy quark 
production in single-tag two-photon events is interesting, because it enables one to compare 
a well-controlled QCD calculation with data at large x and constrain the gluon content of 
the photon at small x. 

Both in the pt and y distributions very little distinguishes the NLO results from the LO. 
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Figure 4: Hadronic and pointlike distributions in pt of the charm quark in NLO for the 
Weizsdcker-Williams (WW) and beamstrahlung cases at a NLC. a) WW case. The upper 
set of histograms is the pointlike compoment, the lower the hadronic. The solid histograms 
represent the LO case, the dashed ones the NLO case, b) Same as in a) but now for the 
(TESLA) beamstrahlung spectrum. 

Assuming a not overly pessimistic charm acceptance, a measurement of the pointhke 
structure function seems feasible at LEP2, and both the pointlike and hadronic structure 
functions are certainly measurable at a future NLC. 
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X 


Events 


range 


range 


lif = Q/2 


/i/ = Q 


11 f = 2Q 


3.2 - 10 


3.2 - 10.0 ■ 10"^ 


3 


2 


2 


1.0-3.2 • 10-3 


25 


21 


20 


3.2 - 10.0 • 10-3 


106 


94 


88 


1.0 - 3.2 • 10^2 


339 


314 


296 


3.2 - 10.0 • 10-2 


919 


884 


857 


1.0 - 3.2 • 10-1 


2350 


2220 


2180 


3.2 - 10.0 ■ 10^1 


1010 


91(3 


815 


10 - 32 


1.0-3.2 • 10-3 


7 


6 


6 


3.2 - 10.0 • 10-3 


63 


58 


55 


1.0 - 3.2 • 10^2 


294 


276 


266 


3.2 - 10.0 • 10-2 


999 


968 


951 


1.0-3.2-10-1 


2910 


2860 


2830 


3.2 - 10.0 • 10-1 


3040 


2880 


2780 


32 - 100 


3.2 - 10.0 • 10-3 


6 


6 


5 


1.0 - 3.2 • 10-2 


54 


51 


50 


3.2 - 10.0 • 10^2 


248 


241 


235 


1.0-3.2-10-1 


825 


814 


808 


3.2-10.0-10-1 


1550 


1520 


1500 


100 - 320 


1.0 - 3.2 ■ 10-2 


5 


5 


5 


3.2 - 10.0 - 10-2 


45 


44 


43 


1.0 - 3.2 - 10-1 


196 


194 


193 


3.2 - 10.0 - 10-1 


511 


509 


508 



Table 1: Number of events with charm quarks due to pointlike photons at LEP2 in ej 
collisions. 
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X 


Events 


range 


range 


fif = Q/2 


/i/ = Q 


/v,/ = 2Q 


3.2 - 10 


3.2 - 10.0 ■ 10"^ 


94 


80 


86 




1.0-3.2 • 10-3 


299 


275 


266 




3.2 - 10.0 • 10-3 


379 


342 


328 




1.0 - 3.2 • 10^2 


328 


290 


271 




3.2 - 10.0 • 10"^ 


192 


159 


144 




1.0 - 3.2 • 10-1 


53 


36 


30 




3.2 - 10.0 ■ 10""^- 


1 








10 - 32 


1.0-3.2 • 10-3 


121 


116 


116 




3.2 - 10.0 • 10-3 


333 


320 


313 




1.0 - 3.2 • 10^2 


410 


385 


372 




3.2 - 10.0 • 10-2 


306 


279 


267 




1.0-3.2-10-1 


111 


95 


88 




3.2 - 10.0 • 10-1 


5 


3 


3 


32 - 100 


3.2 - 10.0 • 10-3 


44 


45 


45 




1.0 - 3.2 • 10-2 


109 


108 


108 




3.2 - 10.0 • 10^2 


112 


110 


110 




1.0-3.2-10-1 


55 


54 


54 




3.2-10.0-10-1 


5 


5 


6 


100 - 320 


1.0 - 3.2 ■ 10-2 


11 


11 


12 




3.2 - 10.0 - 10-2 


23 


24 


24 




1.0 - 3.2 - 10-1 


16 


17 


18 




3.2 - 10.0 - 10-1 


3 


3 


3 



Table 2: Number of events with charm quarks due to hadronic photons at LEP2 in ej 
colhsions. 
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X 


Events 


range 


range 


c(WW) 


5(WW) 


c(Beam) 


6(Beam) 


32 - 100 


1.0-3.2 • 10-=^ 


1850 


120 


30 


2 


3.2 - 10.0 ■ 10-=^ 


2460 


150 


1170 


70 


1.0-3.2-10-2 


2420 


130 


4060 


200 


3.2-10.0- 10-2 


2780 


110 


6510 


260 


1.0-3.2-10-1 


4470 


150 


9590 


320 


3.2-10.0- 10-1 


6740 


60 


11,300 


110 


100 - 320 


1.0-3.2-10-=^ 


1040 


80 








3.2 - 10.0 ■ lO"-'' 


3790 


300 


220 


20 


1.0 - 3.2 - 10-2 


6230 


430 


3150 


200 


3.2-10.0- 10-2 


10,300 


510 


12,200 


580 


1.0-3.2-10^1 


21,700 


830 


34,700 


1320 


3.2-10.0- 10-1 


43,300 


730 


73,800 


1240 


320 - 1000 


1.0-3.2-10-2 


1050 


110 


70 


10 


3.2 - 10.0 ■ 10-- 


2310 


160 


1230 


80 


1.0-3.2 - 10-1 


5610 


260 


6790 


310 


3.2- 10.0- 10-1 


13,300 


400 


21,200 


620 


1000 - 3200 


1.0-3.2-10-2 


70 


10 








3.2-10.0- 10-2 


370 


30 


30 


2 


1.0-3.2-10-1 


1210 


70 


670 


10 


3.2 - 10.0 ■ IQ-i 


3180 


110 


1200 


1(30 



Table 3: Number of events with charm and bottom quarks at a NLC in 67 collisions. 
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Charm 




X 


Events 


Mass 


range 


range 


Uf = 0/2 

H'j ^ / ^ 


lit = O 

H'J ^ 


Uf = 20 


1.3 






2080 


2140 


2230 


1.5 


32 - 100 


1.0 - 3.2 • 10-3 


1880 


1890 


1930 


1.7 






1740 


1750 


1790 


1.3 






346 


356 


370 


1.5 


32 - 100 


3.2 - 10.0 • 10-"^ 


317 


322 


329 


1.7 






296 


305 


313 


1.3 






1060 


1100 


1140 


1.5 


100 - 320 


1.0 - 3.2 • 10-3 


987 


1024 


1058 


1.7 






907 


930 


954 


1.3 






3770 


4080 


4500 


1.5 


32 - 100 


1.0 - 3.2 • 10-3 


3490 


3740 


4090 


1.7 






3240 


3480 


3810 


1.3 






478 


509 


561 


1.5 


32 - 100 


3.2 - 10.0 • 10-"^ 


449 


480 


531 


1.7 






419 


452 


502 


1.3 






829 


916 


1010 


1.5 


100 - 320 


1.0 - 3.2 • 10-3 


774 


854 


943 


1.7 






716 


777 


847 



Table 4: Number of events with charm at small x at a NLC in 67 collisions. We used GRV 
|T5|| MS parton distributions for the top half and ACFGP [|l^ ones for the bottom half. 
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